DISTRIBUTION  STATEMENT  A 

Approved  for  Public  Release 
Distribution  Unlimited 


DRP-1-12 
December  1 993 


Dredging  Research 

T echnical  Notes 


Simulation  of  Wave  Height,  Period,  and  Direction 
Time  Series 

Purpose 

This  technical  note  describes  a  procedure  for  generating  simulated  time 
sequences  of  wave  height,  period,  and  direction  data  corresponding  to  spe¬ 
cific  locations.  The  technique  uses  a  finite-length  wave  record  to  compute 
a  matrix  of  coefficient  multipliers  that  can  be  used  to  generate  arbitrarily 
long  time  sequences  of  simulated  wave  data  which  preserve  the  primary 
statistical  properties  of  the  finite  data  set.  The  procedure  was  developed 
for  simulating  the  20-year  Wave  Information  Study  (WIS)  hindcast  data 
base,  but  can  be  applied  to  any  appropriate  data  sequence. 

Background 

The  goal  of  the  Numerical  Simulation  Techniques  for  Eyaluating  Long- 
Term  Fate  and  Stability'  of  Dredged  Material  Placed  in  Open  Water  work 
unit  is  to  provide  a  simulation  technique  for  determining  how  a  specified 
dredged  material  mound  behaves  over  time.  The  methodology  is  intended 
for  site-designation  investigations  and  is  based  on  long-term  sediment 
transport  simulations  using  local  hydrodynamic  boundary  condition  input 
data.  The  intended  use  of  the  program  is  to  provide  a  systematic  and 
quantifiable  approach  to  analyzing  disposal  site  stability  based  on  local 
environmental  conditions.  If  material  is  eroded  from  a  mound  and  trans¬ 
ported  away  from  the  designated  disposal  site,  the  site  is  classified  as  dis¬ 
persive.  For  locations  predominated  by  strong  wave  and  current  regimes, 
sediment  transport  calculations  based  on  averaged  wave  and  current  data 
may  easily  show  the  site  to  be  dispersive;  however,  if  the  local  environ¬ 
mental  conditions  are  not  severe,  it  may  take  months  or  years  before 
significant  amounts  of  sediment  are  transported  from  the  disposal  site. 

The  ability  to  identify  these  long-term  dispersive  sites  is  especially  impor¬ 
tant  because  the  eroded  material  could  be  transported  into  environmen¬ 
tally  sensitive  areas.  Long-term  dispersion  investigations  require 
knowledge  of  the  local  wave  climate  and  current  field  at  the  site. 
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The  approach  selected  for  disposal  site  analysis  is  a  numerically  based 
sediment  transport  prediction  scheme  using  long-term  local  boundary  con¬ 
ditions  as  input.  These  conditions  represent  those  forcings  that  entrain 
and  transport  sediment,  including  waves,  tides,  and  storms.  In  many 
cases,  these  data  are  not  available  or  are  incomplete  or  too  short  in  dura¬ 
tion  for  long-term  predictions.  To  satisfy  the  ultimate  goals  of  the  work 
unit,  the  intermediate  goal  of  simulating  realistic  boundary  conditions  for 
specific  locations  was  developed.  This  technical  note  addresses  the  genera¬ 
tion  of  arbitrarily  long  time  sequences  of  ocean  wave  data  corresponding 


to  specific  locations.  Subsequent 
components  of  the  current  field. 
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Additional  Information 


Contact  the  authors,  Ms.  Michelle  M.  Thevenot,  (601)  634-3301,  or 
Dr.  Norman  W.  Scheffncr,  (601)  634-3220,  or  the  manager  of  the  Dredging 
Research  Program,  Mr.  E.  Clark  McNair,  (601)  634-2070,  for  additional 
information. 


Introduction 


Unlike  astronomical  tides  that  can  be  determined  with  great  accuracy, 
waves  cannot  be  precisely  predicted.  A  viable  approach  for  generating  re¬ 
alistic  wave  field  boundary  conditions  is  to  develop  a  wave  simulation 
procedure  that  generates  waves  statistically  similar  to  those  known  to 


occur  at  the  site,  that  is,  preserving  seasonality,  directionality,  distribution, 
sequencing,  and  other  characteristics.  A  boundary  condition  database  is 
obtained  by  simulating  an  existing  database  from  which  the  statistical  pa¬ 
rameters  describing  the  intercorrelations  of  wravc  field  parameters  can  be 
computed  (Borgman  and  Scheffner  1991  and  Scheffner  and  Borgman 
1992).  This  procedure  is  somewhat  analogous  to  the  computation  of  har¬ 
monic  constituents  for  a  tidal  record.  The  20-year  WIS  database  (Corson, 
Resio,  and  Vincent  1980  and  Resio,  Vincent,  and  Corson  1982)  w'as  se¬ 
lected  as  the  source  of  wave  parameter  statistics,  from  which  a  matrix  of 
multipliers  is  obtained  that  can  then  be  used  to  generate  time  series  of 
wave  data  which  emulate  the  WIS  hindcasts.  The  advantage  of  the  proce¬ 
dure  is  that  a  finite  length  of  data  can  be  used  to  generate  an  infinitely 
long  data  series  which  is  statistically  similar  to  the  original  n 
liibits  normal  variations  about  its  mean. 
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The  simulation  technique  had  to  meet  two  criteria  for  simulating  time 
sequences  of  wave  data  based  on  finite-length  wave  records.  The  syn¬ 
thetic  record  has  to  be  realistic  with  respect  to  the  statistical  properties  of 

•  i  t 

the  WIS  database,  and  the  methodology  has  to  be  computationally  feasible 
with  respect  to  both  computer  memory  and  computational  speed.  Al¬ 
though  the  intent  of  the  simulation  capability  is  to  provide  times  series 
input  to  the  Long-Term  FATE  model  LTFATE,  the  generated  data  can  be 
used  for  any  application  requiring  wave  data.  The  following  paragraphs 
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describe  the  methodology  required  to  generate  site-specific  wave  field 
time  series. 

Sequence  of  Computations 

In  order  to  generate  arbitrarily  long  time  sequences  of  simulated  wave 
data  which  preserve  the  primary  statistical  properties  of  the  WTS  data, 
two  programs  must  be  run.  First,  the  Height,  Period,  and  Direction  PRE- 
processor  (HPDPRE)  is  run  one  time  to  produce  a  data  summer}'.  The 
data  summary  produced  by  the  preprocessor  can  be  re-used  many  times 
(Figure  1),  with  new  seed  numbers  each  time,  to  obtain  alternate  simula¬ 
tions  with  the  program  Height,  Period,  and  Direction  Simulation 
(HPDSiM).  Both  programs  can  be  run  on  a  top-line  386  personal  com¬ 
puter  with  extended  memory.  This  is  possible  because  both  the  preproces¬ 
sor  HPDPRE  and  the  simulator  HPDSiM  are  optimized  to  work  with 
frequency-domain  methods  and  other  “fast7  algorithms.  Emphasis  has 
been  placed  on  making  the  simulation  program  work  rapidly  because  it 
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Extended  memory  is  not  required  to  run  HPDSIM. 

Example  Application 

The  WIS  station  selected  for  demonstrating  the  capability'  of  the  simula¬ 
tion  package  is  the  Phase’ll  Gulf  of  Mexico  Station  27  (Hubertz  and 


WIS  Data  at  Selected  Site 
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Figure  1.  Flow  diagram  for  creating  multiple  simulations  of  WIS  data 
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Brooks  1989),  located  just  outside  the  entrance  to  Mobile  Bay,  Alabama, 
shown  in  Figure  2.  The  simulated  example  used  is  a  generated  synthetic 
data  series  for  September  1956  through  June  1958,  with  starting  random 
seed  of  67676767  and  is  indexed  as  1957.  Summary  plots  of  the  entire 
year  of  simulated  and  of  W1S  time  series  for  height,  period,  and  direction 
are  shown  in  Figures  3  and  4.  Tables  1  and  2  give  the  computed  values 
of  maximum,  minimum,  average,  and  standard  deviation  for  wave  heights 
and  periods.  The  comparisons  shown  indicate  that  the  simulated  and 
WTS  data  are  similar  with  respect  to  both  magnitudes  and  temporal 
trends,  although  expected  differences  between  the  two  can  be  seen. 

A  more  descriptive  comparison  of  the  simulated  data  to  the  WIS  data 
can  be  seen  through  percent  probability  histograms  that  visually  indicate 
the  percent  of  observations  in  the  data  series  that  fall  within  specified  lim¬ 
its,  for  example,  the  percentage  of  wave  height  observations  between  1.0 
and  1.1  m.  Selected  bands  for  the  histograms  were  0.1-m  increments  for 


Figure  2.  Location  map  for  WIS  Station  27 
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Table  1*  Wave  Height  Comparisons  for  1957 
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Figure  3.  Simulated  wave  height,  period,  and  direction  for  1957 


6 


Technical  Note  DRP-1-12  (December  1993) 


Figure  4.  WIS  wave  height,  period,  and  direction  for  1957 


Technical  Note  DK1M-I2  (December  1993) 


wave  heights,  0.5  sec  for  period,  and  10.0  deg*  for  direction.  These  bands 
can  be  computed  for  both  the  simulated  data  series  and  the  WIS  data 
series  and  overlaid  to  provide  a  visual  comparison  of  both  data  sets.  The 
penod  of  January  is  shown  in  Figure  5  which  illustrates  that  approxi¬ 
mately  30  percent  of  the  January  data  for  both  the  WIS  and  simulated 
time  series  fall  within  the  period  band  of  5.0  to  5.5  sec.  The  plot  provides 
a  comparison  of  the  directionality  of  the  data  that  cannot  be  conveyed  by 
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and  simulated  January  data  series  exhibit  a  bimodal  directionality  with  a 
slight  predominance  in  the  band  from  100  to  120  dev, 


The  variability  of  the  WIS  database  can  be  seen  in  the  direction  histo¬ 
gram  plots  for  June  1956,  1961,  and  1966,  shown  in  Figure  6.  The  year  of 
1956  shows  a  bimodal  distribution  while  1961  and  1966  do  not.  Also,  the 
predominant  direction  of  the  1956  data  is  centered  at  approximately 
90  deg,  but  is  centered  at  nearly  200  deg  in  1966.  A  similar  variability  in 
simulated  data  can  be  demonstrated  by  selecting  different  starting  random 
number  seeds.  Figure  7  presents  three  different  directional  histogram 
plots  for  June  1957,  each  generated  with  a  different  seed.  As  shown,  the 
plots  reflect  the  same  degree  of  variability  as  in  the  WIS  data. 


Summary 


This  statistical  procedure  represents  an  improved  methodology  for  speci¬ 
fying  time  series  boundary  conditions  over  the  use  of  discrete  or  repeated 
historical  data  series.  The  approach  represents  the  stochastic  equivalent  of 
generating  infinitely  long  tidal  series  from  a  finite  number  of  harmonic 
constituents.  The  methodology  gives  the  user  the  ability  to  analyze  the 
variability  of  certain  statistical  parameters  of  an  entire  prototype  database 
such  that  specific  prototype  data  periods  can  be  analyzed  with  respect  to 
their  variability  about  the  mean  conditions.  The  procedure  can  be  used  to 
generate  future  time  series  of  wave  data,  thus  avoiding  the  practice  of 
selecting  design  data  from  past  prototype  conditions  or  of  having  to  re¬ 
peat  fixed  data  series  in  order  to  compose  a  longer  time  series  of  data 
than  is  available.  The  intended  purpose  of  the  database  is  to  generate  infi¬ 
nitely  long,  realistic  wave  data  for  use  as  boundary  conditions  for  predict¬ 
ing  the  long-term  movement  of  dredged  material,  but  the  methodology  is 
recommended  for  any  coastal  application  requiring  either  the  analysis  of 
past  time  series  or  the  generation  of  variable-length  future  time 
boundary  conditions. 


*  To  convert  degrees  to  radians,  multiply  by  0.01745329. 
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Figure  5.  Probability  histogram  comparison  for  January  1957  data 
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